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An  exploratory  investigation  of  measuring  xenon  spectral  emission  lines  in  a  Hall  thruster 
plume  through  optical  imaging  techniques  was  performed  on  the  SPT-100  Hall  thruster.  A 
16  bit  CCD  camera  was  used  with  optical  interference  bandwidth  filters  to  separately  pho¬ 
tograph  neutral  (823  nm)  and  ion  (542  nm)  xenon  emission.  The  resulting  images  allowed 
for  examination  of  the  near-plume  structure  in  terms  of  excited  ions  and  neutral  parti¬ 
cles  for  various  thruster  operating  conditions.  Abel  inversion  analysis  was  applied  to  the 
photographs  to  examine  the  emission  as  a  function  of  radius  with  the  assumption  of  axial 
symmetry  and  an  optically  thin  plasma.  This  experimental  technique  may  provide  results 
useful  for  comparison  with  numerical  near-plume  predictions  that  incorporate  collisional 
radiative  models. 


Nomenclature 

C  constant  that  accounts  for  optical  geometry  parameters  and  physical  constants  in  optical  calibrations 
Id  anode  discharge  current 

fi  fractional  transmittance  of  neutral  density  filter 
TTia  anode  propellant  mass  flow  rate 
rhi)  mass  flow  rate  of  additional  background  gas 
rhc  cathode  propellant  mass  flow  rate 
n  CCD  count  value 

P  vacuum  chamber  pressure 

Q  quantum  efficiency  of  CCD 

r  radial  coordinate  about  thrust  axis 

R  radiance 

Vd  anode  discharge  voltage 

X  axis  in  thruster  coordinate  system  perpendicular  to  the  cathode  plane 

y  axis  in  thruster  coordinate  system  in  the  cathode  plane 

z  axis  in  thruster  coordinate  system  in  thrust  direction 

At  image  exposure  time 

A  wavelength 
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Introduction 


Hall  effect  thruster  plumes  present  a  challenging  environment  for  the  implementation  plasma  diagnostics. 

The  near-plume  region  is  an  especially  difficult  region  to  characterize  where  changes  in  density  and 
temperature  occur  over  short  distances.  Here,  the  physical  presence  of  probes  is  likely  to  perturb  the 
plasma,  affect  thruster  operation,  cause  heating  effects,  and  lead  to  physical  degradation  of  the  probe  itself. 
Even  sizing  a  probe  small  enough  for  good  spatial  resolution  can  be  difficult.  Optical  diagnostics  such  as 
laser-induced  fluorescence^  and  microwave  interferometry^  are  not  encumbered  by  these  limitations  and  have 
proven  useful  for  reliable  near-plume  measurements. 

The  idea  of  using  photography  for  certain  types  of  plasma  characterization  is  intriguing  for  its  poten¬ 
tial  advantages.  Images  provide  data  in  a  virtually  continuous  fashion  over  a  wide  measurement  domain, 
whereas  most  other  diagnostics  typically  measure  at  one  physical  location  at  a  time  and  must  be  physically 
moved  to  create  a  set  of  discrete  data.  The  location  of  the  camera  outside  of  the  chamber  avoids  physical 
disruption  of  the  plasma  and  allows  convenient  access  for  the  operator.  Photographic  images  have  been 
used  to  visually  inspect  basic  changes  in  plume  structure  for  different  thruster  operating  modes. Mod¬ 
ern  high-speed  cameras  can  record  images  on  timescales  that  allow  the  study  of  plasma  oscillation  modes 
in  discharge  channels.^’ ^  Photography  with  optical  filters  has  been  used  to  image  xenon  ion  and  neutral 
emission  separately  in  an  ATON  A53  Hall  thruster.^  This  data  was  transformed  from  line-of-sight  data  into 
cross-sectional  images  of  the  plume  for  two  operational  modes  through  an  inverse  Abel  transformation. 

In  this  study,  a  16  bit  camera  with  optical  interference  bandwidth  filters  was  used  to  examine  the  plasma 
structure  of  the  near-plume  of  the  SPT-100  in  terms  of  optical  emission  of  the  542  nm  xenon  ion  line  and 
the  823  nm  xenon  neutral  line  for  various  operating  conditions.  Optical  emission  is  a  function  of  both 
plasma  density  and  temperature  and  can  be  predicted  by  collisional  radiative  models.  The  goal  of  this  study 
was  to  examine  the  structure  of  the  near-plume  of  the  SPT-100  thruster  at  various  operating  conditions 
and  to  explore  the  use  of  Abel  inversion  analysis  to  deconvolve  the  image  line-of-sight  data  into  a  radially 
dependent  function  with  the  assumption  of  axially  symmetry  and  an  optically  thin  plasma.  The  results  of 
this  experimental  imaging  may  provide  a  useful  comparison  tool  for  numerical  simulations  of  the  near-plume 
region  that  incorporate  collisional  radiative  models. 

Experimental  Apparatus  and  Techniques 


Hall  Effect  Thruster 

A  flight  model  SPT-100  Hall  effect  thruster  (Fig.  1)  was  used  in  this  study.  The  axisymmetric  thruster  is 
equipped  with  two  lanthanum  hexaboride  (LaBe)  cathodes  (only  one  was  used  during  these  tests).  This 
thruster  has  a  conventional  five  magnetic  core  (one  inner,  four  outer)  magnetic  circuit.  Discharge  current 
is  routed  through  the  magnetic  circuit  and  thus  no  extra  power  source  for  the  magnets  is  required.  The 
acceleration  channel  of  the  thruster  has  a  100  mm  outer  diameter,  a  69  mm  inner  diameter,  and  a  channel 
depth  of  25  mm.  For  its  nominal  xenon  operating  condition,  the  thruster  has  been  characterized  to  have  a 
thrust  of  83  mN  with  a  specific  impulse  of  1,600  s,  yielding  an  anode  efficiency  near  50%.^ 

For  this  study,  the  thruster  and  cathode  were  powered  with  commercial  off-the-shelf  Sorenson  power 
supplies  instead  of  the  PPU  used  on-orbit.  A  computer  data  acquisition  system  recorded  the  potential  and 
current  outputs  of  the  power  supplies  used  in  the  thruster  operation  at  a  rate  of  2  Hz.  For  propellant  flow, 
digital  mass  flow  controllers  from  Aera  dispersed  gas  to  the  anode  and  cathode  taking  the  place  of  the 
xenon  flow  control  system  (XFC)  used  on-orbit.  Therefore,  propellant  flow  rate  was  constant  for  each  given 
operating  condition  and  not  controlled  by  discharge  current  feedback  as  with  an  XFC  system. 

Test  Facility 

This  study  utilized  Chamber  1  (Fig.  2(a))  at  the  Air  Force  Research  Laboratory  at  Edwards  Air  Force  Base. 
Chamber  1  is  a  cylindrical  non- magnetic  stainless  steel  vacuum  chamber  2.4  m  in  diameter  and  4.1  m  in 
length.  Pumping  is  provided  by  two  liquid  nitrogen  baffled  (70  K),  1.2  m  flanged  gaseous  helium  two  stage 
cryogenic  (15  K)  vacuum  pumps  with  a  measured  pumping  speed  on  xenon  of  48,500  L/s.  Chamber  pressure 
is  monitored  with  a  hot  filament  ionization  gauge. 
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Figure  1.  The  SPT-100  firing  in  Chamber  1  photographed  by  a  digital  SLR  from  the  same  perspective  as  the 
scientific  images  presented  in  this  study. 


(a)  The  SPT-100  installed  in  Chamber  1  for  the  optical  imaging.  The  (b)  The  SBIG  ST-IOXE  camera  and  optical 
optical  viewport  for  imaging  is  beyond  the  left  photo  border.  components. 


Figure  2.  Experimental  Setup 


Optical  Components 

Camera 

A  Santa  Barbara  Instruments  Group  (SBIG)  ST-IOXE  16  bit  CCD  camera  (Fig.  2(b))was  used  for  this 
imaging  experiment.  It  employs  a  Kodak  KAF-3200E  3  megapixel  CCD  (2184x1472  with  6.8  /am  pixels). 
The  camera  was  designed  for  astronomical  imaging  and  thus  long  exposure  durations  (the  fastest  shutter 
speed  is  0.12  s).  It  interfaces  with  t-mount  optical  components.  The  instrument  has  an  electronic  chiller 
to  reduce  thermal  image  noise,  which  was  set  to  5°C  for  these  measurements.  A  70  mm  focal  length  lens 
was  mounted  to  the  camera  along  with  a  variable  aperture  diaphragm  set  to  a  9.7  mm  diameter  making  for 
an  f- number  of  f/7.2.  Optical  filters  were  stored  and  positioned  with  an  SBIG  CW-8  electronic  filter  wheel. 
The  camera  was  controlled  with  a  laptop  PC.  See  Fig.  4  for  a  schematic  of  the  optical  set  up. 

Filters 

Two  25  mm  diameter  interface  filters  from  Custom  Scientific  were  used  to  isolate  xenon  spectral  emission 
lines.  Singly  charged  ion  emission  was  photographed  through  a  bandwidth  filter  centered  at  542  nm  filter 
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with  a  transmission  distribution  function  having  a  full  width  at  half  maximum  (FWHM)  of  10  nm.  Neutral 
xenon  emission  was  photographed  with  a  filter  centered  at  823  nm  also  with  a  FWHM  of  10  nm.  Besides 
collecting  signal  from  the  823  nm  neutral  line,  this  filter  also  gathers  emission  from  the  828  nm  neutral  line. 
The  high  f- number  of  the  optical  setup  prevented  light  with  high  incidence  angle  from  reaching  the  CCD, 
which  minimized  potential  filter  tuning  effects.  Figure  3  shows  the  filter  transmission  functions  along  with 
the  spectral  radiance  of  the  SPT-100  measured  with  a  spectrometer  through  a  diffuse  optical  collector. 


Figure  3.  Filter  transmission. 


Testing  Methodology 

Optical  Techniques 

The  SPT-100  was  mounted  with  its  centerline  198  cm  from  a  20  cm  dia.  optical  view  port  on  a  lateral  side 
the  vacuum  chamber.  The  camera  was  mounted  on  an  optics  table  outside  the  chamber  and  pointed  through 
the  window  for  a  profile  view  of  the  thruster  and  plume  (Figs.  1  and  4).  A  flocked  graphite  panel  was  placed 
behind  the  thruster  to  serve  as  an  non-reflective  black  background  for  the  images. 

The  shiny  surfaces  of  the  optical  interference  filters  in  conjunction  with  the  other  optical  surfaces  created 
a  mirrored  and  inverted  reflection  of  the  plume  about  the  vertical  centerline  of  the  photo  that  appeared  in  the 
camera  images.  Constraining  the  plume  image  to  one  side  of  the  frame  was  necessary  to  prevent  unwanted 
reflections  from  becoming  superimposed  on  the  primary  plume  image.  Thus,  the  thruster  was  framed  in  the 
camera’s  field  of  view  with  its  plume  confined  to  the  right  half  the  frame  (with  its  downstream  axis  pointing 
rightward).  Unfortunately,  this  solution  to  the  reflection  problem  required  discarding  one  half  of  the  image 
area. 

A  motion  stage  controlled  optical  shade  (Figs.  2(a)  and  4)  was  constructed  from  a  grafoil  (flexible  graphite) 
sheet  covered  with  graphite  felt  and  placed  inside  the  chamber  so  that  certain  regions  of  the  thruster,  which 
varied  drastically  in  brightness,  could  be  selectively  blocked  from  the  view  of  the  camera.  The  shade  could  be 
positioned  so  that  it  blocked  the  cathode  (very  bright  neutral  emission)  or  the  region  of  the  plume  upstream 
of  z=134  mm.  Blocking  bright  regions  of  emission  allowed  for  exposures  of  dimmer  regions  that  more  fully 
utilized  the  bit  depth  of  the  CCD  by  preventing  overexposed  regions  with  pixel  bleeding.  In  some  cases, 
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Figure  4.  Schematic  of  optical  setup 


separate  exposures  for  different  regions  were  combined  to  form  a  composite  photo,  effectively  increasing  the 
dynamic  range  of  the  image.  Three  different  neutral  density  (powers  of  0.75,  2.0,  and  2.5)  were  used  to 
enable  exposure  times  of  between  1  s  and  30  s  and  were  selected  based  on  the  brightness  of  the  region  being 
photographed. 

To  improve  the  signal  to  noise  ratio  in  the  photographs,  10  exposures  were  averaged  for  each  image.  A 
dark  frame  (photo  with  shutter  closed)  for  the  same  exposure  duration  and  same  CCD  temperature  was 
subtracted  to  remove  noise  caused  by  thermal  effects.  Images  were  also  corrected  with  flat  field  images, 
which  are  images  of  a  uniformly  lit  flat  surface.  Images  are  divided  by  the  flat  field  image  to  correct  for 
vignetting  and  variations  in  pixel  response.  Flat  field  corrections  are  important  in  photometric  studies  such 
as  the  present  one  where  vignetting  caused  up  to  a  60%  brightness  variation  between  the  image  center  and 
corner  in  uncorrected  images. 

Calibrations  of  the  optical  system  were  performed  by  photographing  an  integrating  sphere  inside  the 
chamber  illuminated  by  an  Oriel  Model  63355  quartz  tungsten  halogen  lamp  with  a  NIST  traceable  irradiance 
curve.  Calibration  photographs  were  used  to  determine  the  fractional  transmittance  of  each  neutral  density 
filter  at  542  nm  and  823  nm.  The  relative  quantum  efficiency  of  the  CCD  at  542  nm  and  823  nm  was 
also  measured  through  calibration  photos.  With  knowledge  of  these  parameters,  images  taken  with  different 
neutral  density  filters  and  interference  filters  could  be  scaled  relative  to  each  other  in  terms  of  radiance 
according  to  Eq.  1. 


P  Cn 
QXAtfi 

Tests  were  performed  to  measure  pixel  response  with  exposure  time  to  determine  the  point  near  the  upper 
limit  of  pixel  count  range  where  the  CCD  response  became  nonlinear.  All  image  exposure  times  were  selected 
to  ensure  that  pixel  counts  never  exceeded  the  linear  response  limit. 
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Inverse  Abel  Transformation 

These  photographic  images  can  be  considered  to  approximate  a  line-of-sight  data  set  where  plume  emission 
is  integrated  along  the  x-axis  of  the  thruster  coordinate  system.  In  order  to  study  emission  as  a  function 
of  r,  an  Abel  inversion  deconvolution  was  applied  to  image  data.  Besides  the  line-of-sight  assumption,  the 
plasma  was  also  assumed  to  be  optically  thin  and  radially  symmetric  along  the  z-axis  for  a  given  direction 
of  the  ^-axis.  In  this  study,  the  half  of  the  plume  in  the  positive  ^-axis  (non-cathode  side)  was  chosen  for 
Abel  analysis  because  it  was  assumed  to  have  a  higher  degree  of  symmetry  due  to  less  perturbation  from 
cathode. 

The  forward  Abel  transform  for  a  radially  symmetric  function,  /(r)  is  given  by^ 


f{r)rdr 


and  the  analytical  inverse  transformation  is 


(2) 


/w=-- 

TT  Jr  dy  _  ^2 

The  inverse  Abel  transform  is  impractical  to  solve  directly  with  numerical  data  due  to  the  fact  that  it 
requires  differentiation  of  experimental  data,  which  may  be  noisy,  and  contains  a  singularity.^  There  are 
several  types  of  approaches  to  performing  the  inverse  Abel  transform.  Here  the  Fourier  method  described 
in  Ref.  10  was  used  where  the  unknown  radial  function  is  constructed  from  a  series  expansion  of  cosine  terms 
each  multiplied  by  an  amplitude  factor.  The  amplitude  value  for  each  cosine  term  is  determined  through  a 
least  squares  fit  of  the  forward  Abel  transform  equation  with  the  line-of-sight  data.  The  number  of  cosine 
terms  chosen  for  the  expansion  determines  the  amount  of  resolution  and  noise  filtering  in  the  transform.  In 
the  present  study,  50  cosine  terms  were  used. 


Thruster  Operating  Conditions 

Thruster  operating  conditions  for  imaging  were  selected  so  that  the  effects  of  background  pressure,  discharge 
voltage,  and  propellant  flow  rate  on  emission  could  be  characterized.  The  operating  conditions  photographed 
are  summarized  in  Table  1.  To  artificially  increase  chamber  background  pressure,  extra  xenon  gas  was  flowed 
into  the  chamber  through  a  feed  through  at  the  chamber  wall  that  was  5  cm  downstream  from  the  thruster 
exit  plane  and  70  cm  from  the  thruster  center  line.  Discharge  voltage  and  propellant  flow  rates  values  were 
increased  and  decreased  by  30%  of  their  nominal  values.  It  is  important  to  note  that  varying  propellant  flow 
rate  also  affects  this  thruster’s  magnetic  field  because  the  discharge  current  is  routed  through  the  magnet 
coils.  Therefore,  changes  in  the  plume  emission  for  these  cases  are  a  result  of  both  factors. 


Operating  Condition 

Vd  (V) 

Id  (A) 

Power  (W) 

rha  (mg/s) 

rhc  (mg/s) 

rhb  (mg/s) 

P  (Torr-Xe) 

Nominal 

300 

4.61 

1383 

5.15 

0.40 

0.00 

2.28E-05 

High  Back.  Press. 

300 

4.57 

1371 

5.15 

0.40 

4.22 

4.55E-05 

Low  Dis.  Voltage 

210 

4.60 

966 

5.15 

0.40 

0.00 

2.09E-05 

High  Dis.  Voltage 

390 

4.56 

1778 

5.15 

0.40 

0.00 

2.08E-05 

Low  Prop.  Flow  Rate 

300 

2.90 

870 

3.60 

0.28 

0.00 

1.52E-05 

High  Prop.  Flow  Rate 

300 

5.99 

1797 

6.69 

0.51 

0.00 

2.74E-05 

Table  1.  Thruster  operating  specifications. 


Results 

For  each  parameter  variation,  radiance  photograph  sets  are  presented  for  542  nm  and  823  nm  emission 
showing  the  low  and  high  value  parameter  cases  side-by-side  along  with  an  image  showing  the  change  in 
radiance  of  the  high  value  parameter  case  relative  to  the  low  value  case  (Figs.  5,  6,  9,  10,  12,  and  13).  For 
the  delta  images,  gaussian  smoothing  was  applied  to  the  original  images  to  control  noise  amplified  through 
division.  Linear  color  maps  are  applied  to  the  images  to  illustrate  the  variation  in  radiance  near  the  exit 
plane.  However,  radiance  decays  rapidly  beyond  about  a  thruster  diameter  downstream  of  the  exit  plane  so 
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logarithmic  contour  lines  are  also  added  to  delineate  the  shape  and  magnitude  of  the  radiance.  The  lowest 
value  contour  lines  were  sometimes  omitted  where  the  data  was  noisy.  Blocks  were  added  to  the  images  to 
denote  the  channel  geometry.  Only  the  plume  downstream  of  the  exit  plane  was  photographed  so  any  data 
that  appears  to  be  in  the  acceleration  channel  should  be  ignored.  Abel  inversion  transform  analysis  was 
performed  for  the  background  pressure  and  discharge  voltage  variation  cases  at  six  different  z-axis  locations 
(2.5,  5,  10,  20,  40,  and  80  mm).  These  data  are  presented  along  with  the  corresponding  line-of-sight  radiance 
data  (Figs.  7  and  11). 

Figures  5  and  6  show  the  542  and  823  nm  emission,  respectively,  for  the  background  pressure  variation. 
These  photos  show  a  distinct  change  of  the  emission  structure  of  the  plume  for  high  background  pressure 
where  both  ion  and  neutral  emission  decrease  in  the  plume  central  core  beyond  z=b0  mm  and  increase  in  the 
laterally  surrounding  areas.  An  interesting  feature  of  highly  increased  823  nm  emission  (and  542  nm  emission 
to  a  lesser  degree)  is  observed  near  the  exit  plane  outside  of  the  acceleration  channel  on  the  non-cathode  side 
of  the  thruster  {y=7b  mm).  This  feature  is  asymmetric  and  not  observed  on  the  cathode  side  of  the  image. 
With  the  exception  of  the  plume  core,  823  nm  emission  increases  throughout  the  measurement  domain  due  to 
the  higher  density  of  neutral  xenon  atoms  throughout  the  operating  environment.  However,  fig.  5(c)  shows 
a  10%  decrease  in  542  nm  emission  immediately  downstream  of  the  center  magnetic  pole  within  20  mm  of 
the  exit  plane  while  emission  directly  downstream  of  the  acceleration  channel  remains  nearly  constant. 

This  region  of  decreased  emission  in  the  center  of  the  plume  may  indicate  a  decrease  of  interaction  between 
ion  beams  from  opposite  sides  of  the  acceleration  channel  due  to  a  smaller  beamlet  divergence  angle  at  the 
exit  plane.  A  study  of  the  effects  of  background  pressure  on  the  SPT-100  performance  observed  that  thrust 
increases  with  increasing  background  pressure. This  thrust  increase  was  attributed  mostly  to  decreased 
plume  divergence  as  measured  by  probes,  which  may  have  been  caused  by  the  ion  acceleration  zone  shifting 
to  a  more  upstream  location  in  the  discharge  channel  as  observed  in  Ref.  14  for  the  BHT-600  thruster  likely 
due  to  increased  electron  mobility  from  Bohm  diffusion. 

Figure  7(a)  shows  a  similar  ion  emission  radial  profile  (calculated  through  Abel  inversion)  emanating 
from  the  discharge  channel  for  both  background  pressures  through  z=20  mm.  At  40  mm  and  beyond,  the 
lower  emission  at  center  of  the  plume  is  consistent  with  the  observations  made  from  the  line-of-sight  data. 
The  z=2.b  mm  radial  profile  for  the  823  nm  emission  shown  in  Fig.  7(b)  shows  a  distinct  maxima  at  the  edges 
of  the  acceleration  channel  and  a  dip  in  the  center.  Interestingly,  this  distribution  profile  along  the  channel 
exit  is  consistent  with  laser-induced  fluorescence  measurements  of  neutral  density  for  the  SPT-140  described 
in  Ref.  15.  Unlike  the  ion  radial  profiles,  the  neutral  profiles  change  quickly  between  the  two  pressure  cases 
with  downstream  distance.  At  z  =  10  mm  the  local  maxima  from  the  discharge  channel  edges  disappears 
for  the  high  pressure  case  and  the  localized  region  of  high  emission  centered  at  y=7b  mm  and  z=20  mm 
distinctly  appears  as  shown  in  Fig.  6(c). 

The  difference  between  542  and  823  nm  emission  for  the  nominal  operating  condition  is  presented  in 
Fig.  8.  Throughout  the  plume,  the  radiance  of  the  neutral  line  is  stronger.  The  regions  with  the  highest 
relative  542/823  emission  are  the  central  core  of  the  plume  beyond  z=b0  mm  and  near  the  exit  of  the 
acceleration  channel  The  regions  with  the  lowest  relative  emission  are  found  in  the  lateral  areas  outside  of 
the  thruster  plume  where  there  are  very  few  ions. 

Figures  9  and  10  show  the  emission  for  the  discharge  voltage  variation  operating  conditions.  Figure  9(c) 
shows  that  increasing  discharge  voltage  increases  ion  emission  in  the  plume  core  100  mm  downstream  of 
the  exit  plane  while  it  decreases  elsewhere.  This  plume  shape  may  be  indicative  of  the  higher  energy  ions 
undergoing  collisions  from  geometrical  beam  interactions  farther  downstream.  A  similar,  but  more  broad 
collision  region  is  seen  for  the  neutral  emission  in  Fig.  10(c)  which  may  be  from  the  beam  ions  colliding  with 
the  diffuse  neutral  propellant. 

The  radial  ion  emission  profiles  shown  in  Fig.  7(a)  show  a  similar  shape  between  the  two  discharge  voltage 
cases  with  the  low  discharge  voltage  case  emission  significantly  stronger  near  the  exit  plane.  The  magnitude 
of  emission  becomes  nearly  the  same  at  z=80  mm.  A  similar  trend  is  seen  for  the  neutral  emission  in 
Fig.  7(b),  but  the  high  discharge  voltage  emission  in  the  central  plume  is  seen  to  be  greater  for  the  40  and 
80  mm  traces. 

Propellant  flow  rate  effects  (in  combination  with  magnetic  field  effects)  are  shown  in  Figs.  12  and  13. 
Figures  12(b)  and  12(b)  show  that  at  high  propellant  flow  rate,  the  magnitude  of  the  radiance  becomes 
much  greater  and  the  plume  emission  appears  spherical  likely  due  to  a  higher  rate  of  scattering  collisions 
with  the  higher  density  of  unionized  propellant.  This  emission  shape  is  a  drastic  change  from  the  other  test 
conditions  where  the  shape  appears  to  be  influenced  by  geometrical  beam  interactions.  The  total  thruster 
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radiance  increases  dramatically  for  these  high  flow  rate  cases  because  of  the  highly  increased  collisionality 
of  the  plume.  The  increase  in  collisions  happens  mostly  in  the  first  100  mm  of  the  exit  plane. 
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(a)  Nominal  Condition:  542  nm  radiance  (arb.  unit) 
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(b)  Double  Background  Pressure  Condition:  542  nm  radiance 
(arb.  unit) 


150  200mm 


(c)  A  Radiance  (%):  {R542:DP  —  R542:Nom)/ R542:Nom  X  100 

Figure  5.  Ion  emission  radiance  (542  nm)  of  the  nominal  operating  condition  and  the  double  background 
pressure  condition.  Dotted  lines  mark  boundaries  of  composite  photos. 
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150mm  100  50  0  -50  -100  -150 


(a)  Nominal  Conditions:  823  nm  radiance  (arb.  unit) 


(b)  Double  Background  Pressure  Condition:  823  nm  radiance 
(arb.  unit) 

Z 


(c)  A  Radiance  (%):  {R823:DP  —  Rs23-.Nom) / Rs23:Nom  X  100 

Figure  6.  Neutral  emission  radiance  (823  nm)  of  the  nominal  operating  condition  and  double  background 
pressure  condition.  Dotted  lines  mark  boundaries  of  composite  photos. 
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(a)  542  nm  emission 


(b)  823  nm  emission 


Figure  7.  Cross  sections  of  radiance  data  for  the  nominal  and  high  background  pressure  case  at  various  2;  locations  and  their  inverse  Abel  transforms.  Dashed 
lines  represent  acceleration  channel  boundaries. 
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Figure  8.  A  Radiance  (%):  {R^A2:Nom  -  ^823: N om) / ^823 -.Nom  X  100 
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(a)  Vd  =210  V:  542  nm  Radiance  (arb.  unit) 
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Figure  9.  Ion  emission  radiance  (542  nm)  of  the  =210  V  operating  condition  and  the  =390  V  condition. 
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(a)  Vd  =210  V:  823  nm  Radiance  (arb.  unit) 


(b)  Vd  =390  V:  823  nm  Radiance  (arb.  unit) 


(c)  A  Radiance  (%):  {R823:Vd390  —  R823:Vd21o)/R823:Vd210  X 
100 

Figure  10.  Neutral  emission  radiance  (823  nm)  of  the  Vd  =210  V  operating  condition  and  the  Vd  =390  V 
condition.  Missing  image  portion  is  where  optical  shade  blocked  view  of  cathode. 
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Figure  11.  Cross  sections  of  radiance  data  for  the  =210  V  and  =390  V  cases  at  various  2:  locations  and  their  inverse  Abel  transforms.  Dashed  lines 
represent  acceleration  channel  boundaries. 
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Conclusions 


The  imaging  of  a  Hall  thruster  plume  with  spectral  line  filters  shows  distinct  changes  in  the  plume 
emission  structure  from  ions  and  neutral  particles  as  operating  conditions  are  changed.  These  changes 
reflect  variations  in  plume  density  and  temperature.  Photographic  measurements  can  be  used  as  a  convenient 
method  for  qualitatively  examining  plume  structure,  but  may  also  provide  a  source  of  useful  quantitative 
analysis  when  compared  to  collisional  radiative  models  in  numerical  simulations  or  if  they  can  be  correlated 
to  other  diagnostic  measurements  of  density  and/or  temperature. 

The  experimental  techniques  used  in  this  study  showed  good  promise.  Exposures  of  separate  regions  of 
the  plume  with  various  integration  times  and  neutral  density  filters  scaled  remarkably  well  in  magnitude  and 
plume  features  blended  almost  seamlessly  when  combined  for  composite  images.  Overall,  the  Abel  inversion 
of  the  photographic  data  appeared  to  resolve  the  evolution  of  the  plume  from  the  exit  plane  with  good  detail 
showing  effects  that  corresponded  with  the  channel  geometry.  In  a  few  cases,  the  inversion  had  trouble 
resolving  data  near  the  plume  center  and  resulted  in  negative  values,  but  perhaps  using  a  different  range  of 
cosine  frequencies  for  the  solution  or  a  different  inversion  technique  may  alleviate  the  problem. 

In  the  future,  photographic  measurements  may  be  performed  to  match  a  numerically  simulated  data 
set  that  includes  a  collisional  radiative  model  for  model  validation  purposes.  A  more  thorough  optical 
calibration  technique  would  allow  absolute  measurements  of  thruster  radiance,  which  would  add  extra  value 
for  comparison.  Near-held  probe  measurements  of  density  may  be  performed  to  explore  the  possibility  of 
using  photographic  data  to  estimate  density  values  where  temperature  gradients  are  sufficiently  low.  It  may 
be  interesting  to  photograph  these  probes  in  the  near-held  plume  to  examine  the  ehect  their  presence  has 
on  plasma.  The  photographic  techniques  in  this  study  may  also  be  applicable  to  high-speed  photography 
where  time-resolved  measurements  could  be  used  to  study  emission  as  a  function  of  thruster  breathing  mode 
oscillations. 
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Figure  12.  Ion  emission  radiance  (542  nm)  of  the  rh  =-30%  operating  condition  and  the  m  =+30%  condition. 
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Figure  13.  Ion  emission  radiance  (823  nm)  of  the  m  =-30%  operating  condition  and  the  m  =+30%  condition. 
Missing  image  portion  is  where  optical  shade  blocked  view  of  cathode. 
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Introduction 


*  Motivation 

—  Explore  the  possibility  of  using 

photography  as  a  non-intrusive  near¬ 
field  optical  diagnostic 

—  Emission  data  can  be  compared  to 
numerical  simulations  with  collisional 
radiative  models 

*  Summary  of  Experimental  Objectives 

—  Use  a  16  bit  camera  with  optical 
interference  filters  to  image  ion  and 
neutral  emission  lines  (542,  823  nm) 

—  Examine  how  various  operating 
conditions  affect  the  near-field  plume 
emission  of  the  SPT-100 

•  Subtract  photos  to  view  relative 
changes  between  operating 
conditions 

—  Convert  line-of-sight  image  data  to 
radially  dependent  emission  functions 
through  Abel  inversion 
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•  Chamber  1  AFRL  Edwards  AFB 

•  2.4  m  dia.,  4.1  m  length,  stainless  steel 

•  Two  cryogenic  pumps 

-  1 .2  m  dia.,  LN2  baffled  (70  K),  2 

stage  He  (15  K) 

—  SPT-100  nominal  operation  background 

pressure  with  Xe:  2.3x10'®  Torr  (5.55 
mg/s) 


•  Flight  Model  SPT-100 

—  Nominal  Anode  Power:  1350W 
—  Conventional  5  magnetic  core  design 

•  One  inner,  four  outer  connected  in  series 

•  Magnetic  circuit  current  supplied  by  anode 
current 

—  Acceleration  channel:  100  mm  outer  dia.,  69 
mm  inner  dia.,  28  mm  depth 

—  2  lanthanum  hexaboride  (LaBg)  cathodes 
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Optical  Setup 


t 

100  cm 

◄ - ► 


Optical  Shade 


135  mm 

< - ► 


Cathode  1  used 
2  1 


SPT-100 


- ► 

198  cm 
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Optical  Setup  (Continued) 


•  SBIG  ST-1  OXE  Camera 

-  16  bit,  3  MP 

—  f/7.2  optical  setup 

—  Covered  with  shroud  during 
exposures 


Flocked 

Graphite 


*  Chamber  1  Interior  Setup 

—  Thruster  mounted  for  profile  view 

—  Optical  shade  controlled  with  stepper  motor 
driven  stage 

•  Blocks  cathode  and  upstream  plume 

—  Non-reflective  flocked  graphite  panel  for 
background 
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SPT-100  Spectral  Radiance  (arb.  unit) 


Optical  Interference  Filters 


Wavelength  (nm) 
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Optical  Techniques 


Image  Noise  Reduction 

—  Average  10  exposures  for  each 
image 

—  Subtract  dark  frame  from  image 

—  Set  camera’s  electronic  chiller  to  5°  C 

Flat-Field  Correction 

—  Photograph  uniformly  illuminated 
surface 

—  Divide  images  by  flat-field  images  to 
correct  for: 

•  Vignetting  (up  to  60%  brightness 
variation  from  center  to  corner  in 
non-corrected  images) 

•  Non-uniform  pixel  response 

•  Dust  particles  on  CCD 
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Response  (Counts) 


Relative  Radiance  Calibration 


•  Scale  images  through  relative 
radiance  calibration 

—  Integrating  sphere 

—  Quartz  tungsten  lamp  with  NIST 
traceable  irradiance  curve 

•  Calculate  and  Qg23  /  Q542 

•  Characterize  pixel  response  to  ensure 
exposures  are  within  the  linear  regime 


R  = 


=  c 


ens 


/ic^pixel  5'2 


At 


^1  _  ^1  Qx2  ^2  ^^2  fx2 
^2  ^2  Qxi  ^1  fxi 

C  =  combined  physical  and  geometrical  constants 
n  =  CCD  count  value 

=  fractional  transmittance  of  ND  filter  at  a  given  wavelength 
=  CCD  quantum  efficiency  at  a  given  wavelength 
A  =  wavelength 
h  =  Planck  constant 
c  =  speed  of  light 
/tpixei  =  area  of  CCD  pixel 
/\iens  =  area  of  lens 
s'  =  distance  between  lens  and  CCD 
At  =  exposure  time 
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Abel  Inversion  Technique 


Forward 

Transform 


Inverse 
T  ransform 


1  dF(y)  dy 

njy  dy  ^y2  _  j-2 


Inverse  Abel  Transform  is  not  practical 
to  solve  directly  with  numerical  means 

—  Singularity 

—  Requires  derivative  of 
experimental  data 

Use  “Fourier  Method”*  to  approximate 
f(r)  as  a  cosine  expansion 

Nu 

fir)  =  ^  Anfn  (r) 

Ni 

fair)  =  1, 

) 

Choose  frequency  limits  to  balance 
noise  filtering  and  detail  retention 

Least  squares  fit  of  forward  transform 
with  experimental  data  determines 
values  of 


fnir)  =  1-  (-l)"cos  (nn^ 


‘Pretzier,  G.,  Jger,  H.,  Neger,  T,  Philipp,  H.,  and  Woisetschlger,  J.,  “Comparison  of  Different  Methods  of  Abel  Inversion 
Using  Computer  Simulated  and  Experimental  Side-On  Data,"  Zeitschrift  fr  Naturforschung,  Vol.  47a,  1992,  pp.  955-970. 
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Operating  Condition 

V,{V) 

Power 

(W) 

(mg/s) 

m.  (mg/s) 

(mg/s) 

P  (Torr-Xe) 

P  (Pa-Xe) 

Nominal 

300 

4.61 

1383 

5.15 

0.40 

0 

2.28E-05 

3.04E-03 

High  Back.  Press. 

300 

4.57 

1371 

5.15 

0.40 

4.22 

4.55E-05 

6.07E-03 

Low  Dis.  Voltage 

210 

4.60 

966 

5.15 

0.40 

0 

2.09E-05 

2.79E-03 

High  Dis.  Voltage 

390 

4.56 

1778 

5.15 

0.40 

0 

2.08E-05 

2.78E-03 

Low  Prop.  Flow  Rate 

300 

2.90 

870 

3.60 

0.28 

0 

1.52E-05 

2.03E-03 

High  Prop.  Flow  Rate 

300 

5.99 

1797 

6.69 

0.51 

0 

2.74E-05 

3.66E-03 
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150mm  100  50  0  -50  -100  -150 


Comparison  of  Ion  and  Neutral  Emission 

Radiance  (arb.  unit) 


Z 

0  50  100  150  200mm 


o 

lO 


o 

o 


z 

0  50  100  150  200mm 


z 

100  150  200mm 


542  nm  (ion)  823  nm  (neutral)  Afi(%)  =  x 

^823:Wom 

100 

Neutral  emission  is  greater  than  ion  emission  throughout  plume,  especially  at  the 
cathode  (26x  greater  radiance) 

Highest  regions  of  relative  ion  emission  found  at  channel  exits  and  central  plume 
core,  lowest  regions  outside  of  plume  beam 
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150mm  100  50  0  -50  -100  -150 


Z 


542  nm  (Ion)  Emission  Double  BG  Press 

Radiance  (arb.  unit) 


Z 


Z 


50  100  150  200mm 


0  50  100  150  200mm 


Nominal  (P  =  2.3e-5  Torr)  High  Pressure  (P  =  4.6e-5  Torr)  A/?(%)  =  (^542:dp  R542-.Nom)  ^ 

^S42\Nom 

•  Higher  thrust  measured  as  background  pressure  increases  due  mostly  to  less  plume  divergence*  1 00 

•  Asymmetrical  region  of  high  emission  at  z  =  20,  y  =  75  mm 

•  Lower  emission  in  plume  core  and  near  exit  plane,  higher  elsewhere 

—  Maybe  less  beamlet  divergence  from  channel  exits,  from  upstream  shift  of  ion  acceleration  region 

•  Similar  emission  at  channel  exit. 


*Diannant,  K.  D.,  Liang,  R.,  and  Corey,  R.  L.,  “The  Effect  of  Background  Pressure  on  SPT-100  Hall  Thruster  Performance,"  Proceedings  of  the  50th  AIAA/ASME/SAE/ASEE  Joint  Propulsion 
Conference,  2014,  AIAA-201 4-3710. 
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150mm  100  50  0  -50  -100  -150 


Nominal 


Double  BG  Pressure 


A/?(%)  = 


(^823:DP  ^823:JVo?n) 

^823:JVom 


X 


100 

•  Similar,  but  more  intense  asymmetrical  region  of  high  emission  at  z  =  20, 
y  =  75  mm  as  ion  emission 


•  Broader  low  emission  region  in  plume  core  beyond  z  =  50  mm 
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Radiance,  F(y)  (arb.  unit) 


Radial  Profiles:  Double  BG  Press,  vs.  Nominal 


542  nm  (ion)  823  nm  (neutral) 


y,  r(mm) 


y,  r  (mm) 
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Abel  Inversion  of  Radiance.  f(r)  {arb.  unit) 


150mm  100  50  0  -50  -100  -150 


542  nm  (Ion)  Emission  Variations 
Radiance  (arb.  unit) 


z 


z 


50  100  150  200mm 


v^  =  2^oy 


=  390  V 


A/?(%)  = 


(^542:Kd390  ^542:V^d21o) 


^542:Fd210 

100 


X 


•  Greater  emission  in  plume  beyond  z  =  100  mm,  lower  upstream 


Maybe  beam  interaction  is  occurring  farther  downstream  with  the  faster  ions 
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150mm  100  50  0  -50  -100  -150 


Z 

0  50  100  150  200mm 


v^  =  2^oy 


=  390  V 


o 

lO 


A/?(%)  = 


(^823:Kd390  ^823:Kd21o) 

^823:Kd210 

100 


z 

100  150  200mm 


o 

lO 


Similar  changes  compared  to  ion  emission,  but  with  broader  plume  emission  features 
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Radiance,  F(y)  (arb.  unit) 


Radial  Profiles: 


542  nm  (ion) 
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Abel  Inversion  of  Radiance,  f(r)  (arb.  unit) 


150mm  100  50  0  -50  -100  -150 


542  nm  (Ion)  Emission  m  Variations 
Radiance  (arb.  unit) 


z 

0  50  100  150  200mm 


z 

0  50  100  150  200mm 


z 

100  150  200mm 


rh  =  -30%  (m  =  3.88  mg/s) 


rh  =  +30%  {rh  =  7.20  mg/s) 


A/?(%)  = 


(^542:Tn+30%  ^542:Tn-30%) 

^542:m-30% 


X  100 


Dramatic  increase  in  emission  for  high  propellant  flow  rate  (note  color  scale  change) 

—  Highest  increase  directly  downstream  of  center  magnetic  pole  and  outside  of  channel 

—  More  collisions  with  unionized  propellant 
Spherical  emission  shape  different  from  other  cases 

—  Scattering  collisions  diffusing  plume 
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150mm  100  50  0  -50  -100  -150 


823  nm  (Neutral)  Emission  m  Variations 

Radiance  (arb.  unit) 


z 

0  50  100  150  200mm 


rh  =  -30%  (m  =  3.88  mg/s) 


X 


X 


rh  =  +30%  {rh  =  7.20  mg/s) 


A/?(%)  = 


(^823:Tn+30%  ^823:Tn-30%) 

^823:m-30% 


X  100 


z 

0  50  100  150  200mm 


z 

100  150  200mm 


Similar  changes  to  plume  emission  shape  as  seen  for  542  emission 
—  High  increase  in  emission  directly  downstream  of  channel  exit 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  unlimited. 


Conclusions  and  Future  Work 


*  Conclusions 

o  Optical  Technique 

•  High  degree  of  repeatability 

•  Exposures  of  separate  regions  of  plume 
scaled  well  for  composite  image  increasing 
dynamic  range  of  data 

•  Inverse  Abel  inversion  analysis  technique 
performed  well  overall 

—  Resolved  emission  detail  related  to 
acceleration  channel  geometry 

—  Similar  neutral  distribution  shape 
across  channel  as  seen  in  LIF  density 
measurements  SPT-140  near-field 
plume* 

—  A  few  problem  areas  near  r  =  0  where 
f(r)  <  0 

o  Plume  Analysis 

•  Distinct  changes  in  emission  structure  of 
plume  observed  for  various  operating 
conditions 


•  Future  Work 

□  Model  Comparison 

•  Perform  measurements  to  match  a 
numerical  model  of  the  near-field  plume 
that  includes  CRM 

•  Compare  experimental  and  numerical 
results  for  model  validation 

•  Examine  scaling  of  emission  with  density  in 
numerical  results  to  see  if  photographic 
density  estimates  are  appropriate 

□  Probe  Comparison 

•  Directly  measure  density  and/or 
temperature  with  a  probe  in  the  near-field 

•  Compare  results  with  emission  data  to  see 
if  density  can  be  estimated  from  emission 
where  temperature  gradients  are 
adequately  low 

•  Photograph  probe  in  plume  to  examine  the 
perturbation  as  function  of  emission 

□  Apply  experimental  technique  with  high 
speed  camera  to  study  breathing  mode 


*Crofton,  M.  et  al.,  “Neutral  Density  in  the  SPT-140  Near-Field  Plume,"  Proceedings  of  the  33rd  International  Electric 
Propulsion  Conference,  2013,  IEPC-201 3-399. 
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Forward  Abel  Inversion 


hiy)  =  2 


Forward  Abel  Inversion  becomes 


R 

R 

/W  ^ - jd'' 

y 

H(y)  =  2  E  A,  \m—; 

«  =  J  |/j 

dr 


Inverse  Abel  Inversion 


/w-i  j 


Set  each  integral  term  numerically 

R 

r  dr 


r  my)  dy 

Kiy)  = 

A 

Lif) 

J  \/y^-r^ 

m 

} 

\ 

Approximate  f(r)  as  cosine  func.  expansion 

m=  I  AJ„{r) 

F|=iVj 

/oW  =  1  > 


/„(0  =  l-(-ircos|n;r- 


Perform  least  squares  fit  to  fit  amplitude  values  A„ 

n  =  jVj  \  lt=l  /  ic=l 

Vm:  Ni<m< 


N 


Jlr=  1 


‘Pretzier,  G.,  Jger,  H.,  Neger,  T,  Philipp,  H.,  and  Woisetschlger,  J.,  “Comparison  of  Different  Methods  of  Abel  Inversion 
Using  Computer  Simulated  and  Experimental  Side-On  Data,"  Zeitschrift  fr  Naturforschung,  Vol.  47a,  1992,  pp.  955-970. 
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An  exploratory  investigation  of  measuring  xenon  spectral  emission  lines  in  a  Hall  thruster 
plume  through  optical  imaging  techniques  was  performed  on  the  SPT-100  Hall  thruster.  A 
16  bit  CCD  camera  was  used  with  optical  interference  bandwidth  filters  to  separately  pho¬ 
tograph  neutral  (823  nm)  and  ion  (542  nm)  xenon  emission.  The  resulting  images  allowed 
for  examination  of  the  near-plume  structure  in  terms  of  excited  ions  and  neutral  parti¬ 
cles  for  various  thruster  operating  conditions.  Abel  inversion  analysis  was  applied  to  the 
photographs  to  examine  the  emission  as  a  function  of  radius  with  the  assumption  of  axial 
symmetry  and  an  optically  thin  plasma.  This  experimental  technique  may  provide  results 
useful  for  comparison  with  numerical  near-plume  predictions  that  incorporate  collisional 
radiative  models. 


Nomenclature 

C  constant  that  accounts  for  optical  geometry  parameters  and  physical  constants  in  optical  calibrations 
Id  anode  discharge  current 

fi  fractional  transmittance  of  neutral  density  filter 
TTia  anode  propellant  mass  flow  rate 
rhi)  mass  flow  rate  of  additional  background  gas 
rhc  cathode  propellant  mass  flow  rate 
n  CCD  count  value 

P  vacuum  chamber  pressure 

Q  quantum  efficiency  of  CCD 

r  radial  coordinate  about  thrust  axis 

R  radiance 

Vd  anode  discharge  voltage 

X  axis  in  thruster  coordinate  system  perpendicular  to  the  cathode  plane 

y  axis  in  thruster  coordinate  system  in  the  cathode  plane 

z  axis  in  thruster  coordinate  system  in  thrust  direction 

At  image  exposure  time 

A  wavelength 
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Introduction 


Hall  effect  thruster  plumes  present  a  challenging  environment  for  the  implementation  plasma  diagnostics. 

The  near-plume  region  is  an  especially  difficult  region  to  characterize  where  changes  in  density  and 
temperature  occur  over  short  distances.  Here,  the  physical  presence  of  probes  is  likely  to  perturb  the 
plasma,  affect  thruster  operation,  cause  heating  effects,  and  lead  to  physical  degradation  of  the  probe  itself. 
Even  sizing  a  probe  small  enough  for  good  spatial  resolution  can  be  difficult.  Optical  diagnostics  such  as 
laser-induced  fluorescence^  and  microwave  interferometry^  are  not  encumbered  by  these  limitations  and  have 
proven  useful  for  reliable  near-plume  measurements. 

The  idea  of  using  photography  for  certain  types  of  plasma  characterization  is  intriguing  for  its  poten¬ 
tial  advantages.  Images  provide  data  in  a  virtually  continuous  fashion  over  a  wide  measurement  domain, 
whereas  most  other  diagnostics  typically  measure  at  one  physical  location  at  a  time  and  must  be  physically 
moved  to  create  a  set  of  discrete  data.  The  location  of  the  camera  outside  of  the  chamber  avoids  physical 
disruption  of  the  plasma  and  allows  convenient  access  for  the  operator.  Photographic  images  have  been 
used  to  visually  inspect  basic  changes  in  plume  structure  for  different  thruster  operating  modes. Mod¬ 
ern  high-speed  cameras  can  record  images  on  timescales  that  allow  the  study  of  plasma  oscillation  modes 
in  discharge  channels.^’ ^  Photography  with  optical  filters  has  been  used  to  image  xenon  ion  and  neutral 
emission  separately  in  an  ATON  A53  Hall  thruster.^  This  data  was  transformed  from  line-of-sight  data  into 
cross-sectional  images  of  the  plume  for  two  operational  modes  through  an  inverse  Abel  transformation. 

In  this  study,  a  16  bit  camera  with  optical  interference  bandwidth  filters  was  used  to  examine  the  plasma 
structure  of  the  near-plume  of  the  SPT-100  in  terms  of  optical  emission  of  the  542  nm  xenon  ion  line  and 
the  823  nm  xenon  neutral  line  for  various  operating  conditions.  Optical  emission  is  a  function  of  both 
plasma  density  and  temperature  and  can  be  predicted  by  collisional  radiative  models.  The  goal  of  this  study 
was  to  examine  the  structure  of  the  near-plume  of  the  SPT-100  thruster  at  various  operating  conditions 
and  to  explore  the  use  of  Abel  inversion  analysis  to  deconvolve  the  image  line-of-sight  data  into  a  radially 
dependent  function  with  the  assumption  of  axially  symmetry  and  an  optically  thin  plasma.  The  results  of 
this  experimental  imaging  may  provide  a  useful  comparison  tool  for  numerical  simulations  of  the  near-plume 
region  that  incorporate  collisional  radiative  models. 

Experimental  Apparatus  and  Techniques 


Hall  Effect  Thruster 

A  flight  model  SPT-100  Hall  effect  thruster  (Fig.  1)  was  used  in  this  study.  The  axisymmetric  thruster  is 
equipped  with  two  lanthanum  hexaboride  (LaBe)  cathodes  (only  one  was  used  during  these  tests).  This 
thruster  has  a  conventional  five  magnetic  core  (one  inner,  four  outer)  magnetic  circuit.  Discharge  current 
is  routed  through  the  magnetic  circuit  and  thus  no  extra  power  source  for  the  magnets  is  required.  The 
acceleration  channel  of  the  thruster  has  a  100  mm  outer  diameter,  a  69  mm  inner  diameter,  and  a  channel 
depth  of  25  mm.  For  its  nominal  xenon  operating  condition,  the  thruster  has  been  characterized  to  have  a 
thrust  of  83  mN  with  a  specific  impulse  of  1,600  s,  yielding  an  anode  efficiency  near  50%.^ 

For  this  study,  the  thruster  and  cathode  were  powered  with  commercial  off-the-shelf  Sorenson  power 
supplies  instead  of  the  PPU  used  on-orbit.  A  computer  data  acquisition  system  recorded  the  potential  and 
current  outputs  of  the  power  supplies  used  in  the  thruster  operation  at  a  rate  of  2  Hz.  For  propellant  flow, 
digital  mass  flow  controllers  from  Aera  dispersed  gas  to  the  anode  and  cathode  taking  the  place  of  the 
xenon  flow  control  system  (XFC)  used  on-orbit.  Therefore,  propellant  flow  rate  was  constant  for  each  given 
operating  condition  and  not  controlled  by  discharge  current  feedback  as  with  an  XFC  system. 

Test  Facility 

This  study  utilized  Chamber  1  (Fig.  2(a))  at  the  Air  Force  Research  Laboratory  at  Edwards  Air  Force  Base. 
Chamber  1  is  a  cylindrical  non- magnetic  stainless  steel  vacuum  chamber  2.4  m  in  diameter  and  4.1  m  in 
length.  Pumping  is  provided  by  two  liquid  nitrogen  baffled  (70  K),  1.2  m  flanged  gaseous  helium  two  stage 
cryogenic  (15  K)  vacuum  pumps  with  a  measured  pumping  speed  on  xenon  of  48,500  L/s.  Chamber  pressure 
is  monitored  with  a  hot  filament  ionization  gauge. 
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Figure  1.  The  SPT-100  firing  in  Chamber  1  photographed  by  a  digital  SLR  from  the  same  perspective  as  the 
scientific  images  presented  in  this  study. 


(a)  The  SPT-100  installed  in  Chamber  1  for  the  optical  imaging.  The  (b)  The  SBIG  ST-IOXE  camera  and  optical 
optical  viewport  for  imaging  is  beyond  the  left  photo  border.  components. 


Figure  2.  Experimental  Setup 


Optical  Components 

Camera 

A  Santa  Barbara  Instruments  Group  (SBIG)  ST-IOXE  16  bit  CCD  camera  (Fig.  2(b))was  used  for  this 
imaging  experiment.  It  employs  a  Kodak  KAF-3200E  3  megapixel  CCD  (2184x1472  with  6.8  /am  pixels). 
The  camera  was  designed  for  astronomical  imaging  and  thus  long  exposure  durations  (the  fastest  shutter 
speed  is  0.12  s).  It  interfaces  with  t-mount  optical  components.  The  instrument  has  an  electronic  chiller 
to  reduce  thermal  image  noise,  which  was  set  to  5°C  for  these  measurements.  A  70  mm  focal  length  lens 
was  mounted  to  the  camera  along  with  a  variable  aperture  diaphragm  set  to  a  9.7  mm  diameter  making  for 
an  f- number  of  f/7.2.  Optical  filters  were  stored  and  positioned  with  an  SBIG  CW-8  electronic  filter  wheel. 
The  camera  was  controlled  with  a  laptop  PC.  See  Fig.  4  for  a  schematic  of  the  optical  set  up. 

Filters 

Two  25  mm  diameter  interface  filters  from  Custom  Scientific  were  used  to  isolate  xenon  spectral  emission 
lines.  Singly  charged  ion  emission  was  photographed  through  a  bandwidth  filter  centered  at  542  nm  filter 
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with  a  transmission  distribution  function  having  a  full  width  at  half  maximum  (FWHM)  of  10  nm.  Neutral 
xenon  emission  was  photographed  with  a  filter  centered  at  823  nm  also  with  a  FWHM  of  10  nm.  Besides 
collecting  signal  from  the  823  nm  neutral  line,  this  filter  also  gathers  emission  from  the  828  nm  neutral  line. 
The  high  f- number  of  the  optical  setup  prevented  light  with  high  incidence  angle  from  reaching  the  CCD, 
which  minimized  potential  filter  tuning  effects.  Figure  3  shows  the  filter  transmission  functions  along  with 
the  spectral  radiance  of  the  SPT-100  measured  with  a  spectrometer  through  a  diffuse  optical  collector. 


Figure  3.  Filter  transmission. 


Testing  Methodology 

Optical  Techniques 

The  SPT-100  was  mounted  with  its  centerline  198  cm  from  a  20  cm  dia.  optical  view  port  on  a  lateral  side 
the  vacuum  chamber.  The  camera  was  mounted  on  an  optics  table  outside  the  chamber  and  pointed  through 
the  window  for  a  profile  view  of  the  thruster  and  plume  (Figs.  1  and  4).  A  flocked  graphite  panel  was  placed 
behind  the  thruster  to  serve  as  an  non-reflective  black  background  for  the  images. 

The  shiny  surfaces  of  the  optical  interference  filters  in  conjunction  with  the  other  optical  surfaces  created 
a  mirrored  and  inverted  reflection  of  the  plume  about  the  vertical  centerline  of  the  photo  that  appeared  in  the 
camera  images.  Constraining  the  plume  image  to  one  side  of  the  frame  was  necessary  to  prevent  unwanted 
reflections  from  becoming  superimposed  on  the  primary  plume  image.  Thus,  the  thruster  was  framed  in  the 
camera’s  field  of  view  with  its  plume  confined  to  the  right  half  the  frame  (with  its  downstream  axis  pointing 
rightward).  Unfortunately,  this  solution  to  the  reflection  problem  required  discarding  one  half  of  the  image 
area. 

A  motion  stage  controlled  optical  shade  (Figs.  2(a)  and  4)  was  constructed  from  a  grafoil  (flexible  graphite) 
sheet  covered  with  graphite  felt  and  placed  inside  the  chamber  so  that  certain  regions  of  the  thruster,  which 
varied  drastically  in  brightness,  could  be  selectively  blocked  from  the  view  of  the  camera.  The  shade  could  be 
positioned  so  that  it  blocked  the  cathode  (very  bright  neutral  emission)  or  the  region  of  the  plume  upstream 
of  z=134  mm.  Blocking  bright  regions  of  emission  allowed  for  exposures  of  dimmer  regions  that  more  fully 
utilized  the  bit  depth  of  the  CCD  by  preventing  overexposed  regions  with  pixel  bleeding.  In  some  cases, 
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Figure  4.  Schematic  of  optical  setup 


separate  exposures  for  different  regions  were  combined  to  form  a  composite  photo,  effectively  increasing  the 
dynamic  range  of  the  image.  Three  different  neutral  density  (powers  of  0.75,  2.0,  and  2.5)  were  used  to 
enable  exposure  times  of  between  1  s  and  30  s  and  were  selected  based  on  the  brightness  of  the  region  being 
photographed. 

To  improve  the  signal  to  noise  ratio  in  the  photographs,  10  exposures  were  averaged  for  each  image.  A 
dark  frame  (photo  with  shutter  closed)  for  the  same  exposure  duration  and  same  CCD  temperature  was 
subtracted  to  remove  noise  caused  by  thermal  effects.  Images  were  also  corrected  with  flat  field  images, 
which  are  images  of  a  uniformly  lit  flat  surface.  Images  are  divided  by  the  flat  field  image  to  correct  for 
vignetting  and  variations  in  pixel  response.  Flat  field  corrections  are  important  in  photometric  studies  such 
as  the  present  one  where  vignetting  caused  up  to  a  60%  brightness  variation  between  the  image  center  and 
corner  in  uncorrected  images. 

Calibrations  of  the  optical  system  were  performed  by  photographing  an  integrating  sphere  inside  the 
chamber  illuminated  by  an  Oriel  Model  63355  quartz  tungsten  halogen  lamp  with  a  NIST  traceable  irradiance 
curve.  Calibration  photographs  were  used  to  determine  the  fractional  transmittance  of  each  neutral  density 
filter  at  542  nm  and  823  nm.  The  relative  quantum  efficiency  of  the  CCD  at  542  nm  and  823  nm  was 
also  measured  through  calibration  photos.  With  knowledge  of  these  parameters,  images  taken  with  different 
neutral  density  filters  and  interference  filters  could  be  scaled  relative  to  each  other  in  terms  of  radiance 
according  to  Eq.  1. 


P  Cn 
QXAtfi 

Tests  were  performed  to  measure  pixel  response  with  exposure  time  to  determine  the  point  near  the  upper 
limit  of  pixel  count  range  where  the  CCD  response  became  nonlinear.  All  image  exposure  times  were  selected 
to  ensure  that  pixel  counts  never  exceeded  the  linear  response  limit. 
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Inverse  Abel  Transformation 

These  photographic  images  can  be  considered  to  approximate  a  line-of-sight  data  set  where  plume  emission 
is  integrated  along  the  x-axis  of  the  thruster  coordinate  system.  In  order  to  study  emission  as  a  function 
of  r,  an  Abel  inversion  deconvolution  was  applied  to  image  data.  Besides  the  line-of-sight  assumption,  the 
plasma  was  also  assumed  to  be  optically  thin  and  radially  symmetric  along  the  z-axis  for  a  given  direction 
of  the  ^-axis.  In  this  study,  the  half  of  the  plume  in  the  positive  ^-axis  (non-cathode  side)  was  chosen  for 
Abel  analysis  because  it  was  assumed  to  have  a  higher  degree  of  symmetry  due  to  less  perturbation  from 
cathode. 

The  forward  Abel  transform  for  a  radially  symmetric  function,  /(r)  is  given  by^ 


f{r)rdr 


and  the  analytical  inverse  transformation  is 


(2) 


/w=-- 

TT  Jr  dy  _  ^2 

The  inverse  Abel  transform  is  impractical  to  solve  directly  with  numerical  data  due  to  the  fact  that  it 
requires  differentiation  of  experimental  data,  which  may  be  noisy,  and  contains  a  singularity.^  There  are 
several  types  of  approaches  to  performing  the  inverse  Abel  transform.  Here  the  Fourier  method  described 
in  Ref.  10  was  used  where  the  unknown  radial  function  is  constructed  from  a  series  expansion  of  cosine  terms 
each  multiplied  by  an  amplitude  factor.  The  amplitude  value  for  each  cosine  term  is  determined  through  a 
least  squares  fit  of  the  forward  Abel  transform  equation  with  the  line-of-sight  data.  The  number  of  cosine 
terms  chosen  for  the  expansion  determines  the  amount  of  resolution  and  noise  filtering  in  the  transform.  In 
the  present  study,  50  cosine  terms  were  used. 


Thruster  Operating  Conditions 

Thruster  operating  conditions  for  imaging  were  selected  so  that  the  effects  of  background  pressure,  discharge 
voltage,  and  propellant  flow  rate  on  emission  could  be  characterized.  The  operating  conditions  photographed 
are  summarized  in  Table  1.  To  artificially  increase  chamber  background  pressure,  extra  xenon  gas  was  flowed 
into  the  chamber  through  a  feed  through  at  the  chamber  wall  that  was  5  cm  downstream  from  the  thruster 
exit  plane  and  70  cm  from  the  thruster  center  line.  Discharge  voltage  and  propellant  flow  rates  values  were 
increased  and  decreased  by  30%  of  their  nominal  values.  It  is  important  to  note  that  varying  propellant  flow 
rate  also  affects  this  thruster’s  magnetic  field  because  the  discharge  current  is  routed  through  the  magnet 
coils.  Therefore,  changes  in  the  plume  emission  for  these  cases  are  a  result  of  both  factors. 


Operating  Condition 

Vd  (V) 

Id  (A) 

Power  (W) 

rha  (mg/s) 

rhc  (mg/s) 

rhb  (mg/s) 

P  (Torr-Xe) 

Nominal 

300 

4.61 

1383 

5.15 

0.40 

0.00 

2.28E-05 

High  Back.  Press. 

300 

4.57 

1371 

5.15 

0.40 

4.22 

4.55E-05 

Low  Dis.  Voltage 

210 

4.60 

966 

5.15 

0.40 

0.00 

2.09E-05 

High  Dis.  Voltage 

390 

4.56 

1778 

5.15 

0.40 

0.00 

2.08E-05 

Low  Prop.  Flow  Rate 

300 

2.90 

870 

3.60 

0.28 

0.00 

1.52E-05 

High  Prop.  Flow  Rate 

300 

5.99 

1797 

6.69 

0.51 

0.00 

2.74E-05 

Table  1.  Thruster  operating  specifications. 


Results 

For  each  parameter  variation,  radiance  photograph  sets  are  presented  for  542  nm  and  823  nm  emission 
showing  the  low  and  high  value  parameter  cases  side-by-side  along  with  an  image  showing  the  change  in 
radiance  of  the  high  value  parameter  case  relative  to  the  low  value  case  (Figs.  5,  6,  9,  10,  12,  and  13).  For 
the  delta  images,  gaussian  smoothing  was  applied  to  the  original  images  to  control  noise  amplified  through 
division.  Linear  color  maps  are  applied  to  the  images  to  illustrate  the  variation  in  radiance  near  the  exit 
plane.  However,  radiance  decays  rapidly  beyond  about  a  thruster  diameter  downstream  of  the  exit  plane  so 
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logarithmic  contour  lines  are  also  added  to  delineate  the  shape  and  magnitude  of  the  radiance.  The  lowest 
value  contour  lines  were  sometimes  omitted  where  the  data  was  noisy.  Blocks  were  added  to  the  images  to 
denote  the  channel  geometry.  Only  the  plume  downstream  of  the  exit  plane  was  photographed  so  any  data 
that  appears  to  be  in  the  acceleration  channel  should  be  ignored.  Abel  inversion  transform  analysis  was 
performed  for  the  background  pressure  and  discharge  voltage  variation  cases  at  six  different  z-axis  locations 
(2.5,  5,  10,  20,  40,  and  80  mm).  These  data  are  presented  along  with  the  corresponding  line-of-sight  radiance 
data  (Figs.  7  and  11). 

Figures  5  and  6  show  the  542  and  823  nm  emission,  respectively,  for  the  background  pressure  variation. 
These  photos  show  a  distinct  change  of  the  emission  structure  of  the  plume  for  high  background  pressure 
where  both  ion  and  neutral  emission  decrease  in  the  plume  central  core  beyond  z=b0  mm  and  increase  in  the 
laterally  surrounding  areas.  An  interesting  feature  of  highly  increased  823  nm  emission  (and  542  nm  emission 
to  a  lesser  degree)  is  observed  near  the  exit  plane  outside  of  the  acceleration  channel  on  the  non-cathode  side 
of  the  thruster  {y=7b  mm).  This  feature  is  asymmetric  and  not  observed  on  the  cathode  side  of  the  image. 
With  the  exception  of  the  plume  core,  823  nm  emission  increases  throughout  the  measurement  domain  due  to 
the  higher  density  of  neutral  xenon  atoms  throughout  the  operating  environment.  However,  fig.  5(c)  shows 
a  10%  decrease  in  542  nm  emission  immediately  downstream  of  the  center  magnetic  pole  within  20  mm  of 
the  exit  plane  while  emission  directly  downstream  of  the  acceleration  channel  remains  nearly  constant. 

This  region  of  decreased  emission  in  the  center  of  the  plume  may  indicate  a  decrease  of  interaction  between 
ion  beams  from  opposite  sides  of  the  acceleration  channel  due  to  a  smaller  beamlet  divergence  angle  at  the 
exit  plane.  A  study  of  the  effects  of  background  pressure  on  the  SPT-100  performance  observed  that  thrust 
increases  with  increasing  background  pressure. This  thrust  increase  was  attributed  mostly  to  decreased 
plume  divergence  as  measured  by  probes,  which  may  have  been  caused  by  the  ion  acceleration  zone  shifting 
to  a  more  upstream  location  in  the  discharge  channel  as  observed  in  Ref.  14  for  the  BHT-600  thruster  likely 
due  to  increased  electron  mobility  from  Bohm  diffusion. 

Figure  7(a)  shows  a  similar  ion  emission  radial  profile  (calculated  through  Abel  inversion)  emanating 
from  the  discharge  channel  for  both  background  pressures  through  z=20  mm.  At  40  mm  and  beyond,  the 
lower  emission  at  center  of  the  plume  is  consistent  with  the  observations  made  from  the  line-of-sight  data. 
The  z=2.b  mm  radial  profile  for  the  823  nm  emission  shown  in  Fig.  7(b)  shows  a  distinct  maxima  at  the  edges 
of  the  acceleration  channel  and  a  dip  in  the  center.  Interestingly,  this  distribution  profile  along  the  channel 
exit  is  consistent  with  laser-induced  fluorescence  measurements  of  neutral  density  for  the  SPT-140  described 
in  Ref.  15.  Unlike  the  ion  radial  profiles,  the  neutral  profiles  change  quickly  between  the  two  pressure  cases 
with  downstream  distance.  At  z  =  10  mm  the  local  maxima  from  the  discharge  channel  edges  disappears 
for  the  high  pressure  case  and  the  localized  region  of  high  emission  centered  at  y=7b  mm  and  z=20  mm 
distinctly  appears  as  shown  in  Fig.  6(c). 

The  difference  between  542  and  823  nm  emission  for  the  nominal  operating  condition  is  presented  in 
Fig.  8.  Throughout  the  plume,  the  radiance  of  the  neutral  line  is  stronger.  The  regions  with  the  highest 
relative  542/823  emission  are  the  central  core  of  the  plume  beyond  z=b0  mm  and  near  the  exit  of  the 
acceleration  channel  The  regions  with  the  lowest  relative  emission  are  found  in  the  lateral  areas  outside  of 
the  thruster  plume  where  there  are  very  few  ions. 

Figures  9  and  10  show  the  emission  for  the  discharge  voltage  variation  operating  conditions.  Figure  9(c) 
shows  that  increasing  discharge  voltage  increases  ion  emission  in  the  plume  core  100  mm  downstream  of 
the  exit  plane  while  it  decreases  elsewhere.  This  plume  shape  may  be  indicative  of  the  higher  energy  ions 
undergoing  collisions  from  geometrical  beam  interactions  farther  downstream.  A  similar,  but  more  broad 
collision  region  is  seen  for  the  neutral  emission  in  Fig.  10(c)  which  may  be  from  the  beam  ions  colliding  with 
the  diffuse  neutral  propellant. 

The  radial  ion  emission  profiles  shown  in  Fig.  7(a)  show  a  similar  shape  between  the  two  discharge  voltage 
cases  with  the  low  discharge  voltage  case  emission  significantly  stronger  near  the  exit  plane.  The  magnitude 
of  emission  becomes  nearly  the  same  at  z=80  mm.  A  similar  trend  is  seen  for  the  neutral  emission  in 
Fig.  7(b),  but  the  high  discharge  voltage  emission  in  the  central  plume  is  seen  to  be  greater  for  the  40  and 
80  mm  traces. 

Propellant  flow  rate  effects  (in  combination  with  magnetic  field  effects)  are  shown  in  Figs.  12  and  13. 
Figures  12(b)  and  12(b)  show  that  at  high  propellant  flow  rate,  the  magnitude  of  the  radiance  becomes 
much  greater  and  the  plume  emission  appears  spherical  likely  due  to  a  higher  rate  of  scattering  collisions 
with  the  higher  density  of  unionized  propellant.  This  emission  shape  is  a  drastic  change  from  the  other  test 
conditions  where  the  shape  appears  to  be  influenced  by  geometrical  beam  interactions.  The  total  thruster 
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radiance  increases  dramatically  for  these  high  flow  rate  cases  because  of  the  highly  increased  collisionality 
of  the  plume.  The  increase  in  collisions  happens  mostly  in  the  first  100  mm  of  the  exit  plane. 
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(a)  Nominal  Condition:  542  nm  radiance  (arb.  unit) 
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(b)  Double  Background  Pressure  Condition:  542  nm  radiance 
(arb.  unit) 
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(c)  A  Radiance  (%):  {R542:DP  —  R542:Nom)/ R542:Nom  X  100 

Figure  5.  Ion  emission  radiance  (542  nm)  of  the  nominal  operating  condition  and  the  double  background 
pressure  condition.  Dotted  lines  mark  boundaries  of  composite  photos. 
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(a)  Nominal  Conditions:  823  nm  radiance  (arb.  unit) 


(b)  Double  Background  Pressure  Condition:  823  nm  radiance 
(arb.  unit) 
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(c)  A  Radiance  (%):  {R823:DP  —  Rs23-.Nom) / Rs23:Nom  X  100 

Figure  6.  Neutral  emission  radiance  (823  nm)  of  the  nominal  operating  condition  and  double  background 
pressure  condition.  Dotted  lines  mark  boundaries  of  composite  photos. 
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(a)  542  nm  emission 


(b)  823  nm  emission 


Figure  7.  Cross  sections  of  radiance  data  for  the  nominal  and  high  background  pressure  case  at  various  2;  locations  and  their  inverse  Abel  transforms.  Dashed 
lines  represent  acceleration  channel  boundaries. 


Abel  Inversion  of  Radiance,  f(r}  (arb.  unit) 


z 


Figure  8.  A  Radiance  (%):  {R^A2:Nom  -  ^823: N om) / ^823 -.Nom  X  100 
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150mm  100  50  0  -50  -100  -150 


(a)  Vd  =210  V:  542  nm  Radiance  (arb.  unit) 


(b)  Vd  =390  V:  542  nm  Radiance  (arb.  unit) 


(c)  A  Radiance  (%):  {R542:Vd390  —  R542:Vd21o)/R542:Vd210  X 
100 


Figure  9.  Ion  emission  radiance  (542  nm)  of  the  =210  V  operating  condition  and  the  =390  V  condition. 
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150mm  100  50  0  -50  -100  -150 


z 


z 


(a)  Vd  =210  V:  823  nm  Radiance  (arb.  unit) 


(b)  Vd  =390  V:  823  nm  Radiance  (arb.  unit) 


(c)  A  Radiance  (%):  {R823:Vd390  —  R823:Vd21o)/R823:Vd210  X 
100 

Figure  10.  Neutral  emission  radiance  (823  nm)  of  the  Vd  =210  V  operating  condition  and  the  Vd  =390  V 
condition.  Missing  image  portion  is  where  optical  shade  blocked  view  of  cathode. 
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Figure  11.  Cross  sections  of  radiance  data  for  the  =210  V  and  =390  V  cases  at  various  2:  locations  and  their  inverse  Abel  transforms.  Dashed  lines 
represent  acceleration  channel  boundaries. 
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Conclusions 


The  imaging  of  a  Hall  thruster  plume  with  spectral  line  filters  shows  distinct  changes  in  the  plume 
emission  structure  from  ions  and  neutral  particles  as  operating  conditions  are  changed.  These  changes 
reflect  variations  in  plume  density  and  temperature.  Photographic  measurements  can  be  used  as  a  convenient 
method  for  qualitatively  examining  plume  structure,  but  may  also  provide  a  source  of  useful  quantitative 
analysis  when  compared  to  collisional  radiative  models  in  numerical  simulations  or  if  they  can  be  correlated 
to  other  diagnostic  measurements  of  density  and/or  temperature. 

The  experimental  techniques  used  in  this  study  showed  good  promise.  Exposures  of  separate  regions  of 
the  plume  with  various  integration  times  and  neutral  density  filters  scaled  remarkably  well  in  magnitude  and 
plume  features  blended  almost  seamlessly  when  combined  for  composite  images.  Overall,  the  Abel  inversion 
of  the  photographic  data  appeared  to  resolve  the  evolution  of  the  plume  from  the  exit  plane  with  good  detail 
showing  effects  that  corresponded  with  the  channel  geometry.  In  a  few  cases,  the  inversion  had  trouble 
resolving  data  near  the  plume  center  and  resulted  in  negative  values,  but  perhaps  using  a  different  range  of 
cosine  frequencies  for  the  solution  or  a  different  inversion  technique  may  alleviate  the  problem. 

In  the  future,  photographic  measurements  may  be  performed  to  match  a  numerically  simulated  data 
set  that  includes  a  collisional  radiative  model  for  model  validation  purposes.  A  more  thorough  optical 
calibration  technique  would  allow  absolute  measurements  of  thruster  radiance,  which  would  add  extra  value 
for  comparison.  Near-held  probe  measurements  of  density  may  be  performed  to  explore  the  possibility  of 
using  photographic  data  to  estimate  density  values  where  temperature  gradients  are  sufficiently  low.  It  may 
be  interesting  to  photograph  these  probes  in  the  near-held  plume  to  examine  the  ehect  their  presence  has 
on  plasma.  The  photographic  techniques  in  this  study  may  also  be  applicable  to  high-speed  photography 
where  time-resolved  measurements  could  be  used  to  study  emission  as  a  function  of  thruster  breathing  mode 
oscillations. 
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150mm  100  50  0  -50  -100  -150 


z 


z 


(a)  m  =-30%:  542  nm  Radiance  (arb.  unit) 


(b)  rh  =+30%:  542  nm  Radiance  (arb.  unit)  Note:  Color  Scale 
Change 

Z 


(c)  A  Radiance  (%):  (^542:m+30% 

^542:m-30%)/ ^542  :m  — 30%  ^  100 

Figure  12.  Ion  emission  radiance  (542  nm)  of  the  rh  =-30%  operating  condition  and  the  m  =+30%  condition. 
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150mm  100  50  0  -50  -100  -150 


(a)  m 


=-30%:  823  nm  Radiance  (arb.  unit) 


(b)  m  =+30%:  823  nm  Radiance  (arb.  unit)  Note:  Color  Scale 
Change 

Z 


(c)  A  Radiance  (%):  (^823:m+30% 

^823:m-30%)/^823:m-30%  X  100 

Figure  13.  Ion  emission  radiance  (823  nm)  of  the  m  =-30%  operating  condition  and  the  m  =+30%  condition. 
Missing  image  portion  is  where  optical  shade  blocked  view  of  cathode. 
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Introduction 


*  Motivation 

—  Explore  the  possibility  of  using 

photography  as  a  non-intrusive  near¬ 
field  optical  diagnostic 

—  Emission  data  can  be  compared  to 
numerical  simulations  with  collisional 
radiative  models 

*  Summary  of  Experimental  Objectives 

—  Use  a  16  bit  camera  with  optical 
interference  filters  to  image  ion  and 
neutral  emission  lines  (542,  823  nm) 

—  Examine  how  various  operating 
conditions  affect  the  near-field  plume 
emission  of  the  SPT-100 

•  Subtract  photos  to  view  relative 
changes  between  operating 
conditions 

—  Convert  line-of-sight  image  data  to 
radially  dependent  emission  functions 
through  Abel  inversion 
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•  Chamber  1  AFRL  Edwards  AFB 

•  2.4  m  dia.,  4.1  m  length,  stainless  steel 

•  Two  cryogenic  pumps 

-  1 .2  m  dia.,  LN2  baffled  (70  K),  2 

stage  He  (15  K) 

—  SPT-100  nominal  operation  background 

pressure  with  Xe:  2.3x10'®  Torr  (5.55 
mg/s) 


•  Flight  Model  SPT-100 

—  Nominal  Anode  Power:  1350W 
—  Conventional  5  magnetic  core  design 

•  One  inner,  four  outer  connected  in  series 

•  Magnetic  circuit  current  supplied  by  anode 
current 

—  Acceleration  channel:  100  mm  outer  dia.,  69 
mm  inner  dia.,  28  mm  depth 

—  2  lanthanum  hexaboride  (LaBg)  cathodes 
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Optical  Setup 


t 

100  cm 

◄ - ► 


Optical  Shade 


135  mm 

< - ► 


Cathode  1  used 
2  1 


SPT-100 


- ► 

198  cm 
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Optical  Setup  (Continued) 


•  SBIG  ST-1  OXE  Camera 

-  16  bit,  3  MP 

—  f/7.2  optical  setup 

—  Covered  with  shroud  during 
exposures 


Flocked 

Graphite 


*  Chamber  1  Interior  Setup 

—  Thruster  mounted  for  profile  view 

—  Optical  shade  controlled  with  stepper  motor 
driven  stage 

•  Blocks  cathode  and  upstream  plume 

—  Non-reflective  flocked  graphite  panel  for 
background 
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SPT-100  Spectral  Radiance  (arb.  unit) 


Optical  Interference  Filters 


Wavelength  (nm) 
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Optical  Techniques 


Image  Noise  Reduction 

—  Average  10  exposures  for  each 
image 

—  Subtract  dark  frame  from  image 

—  Set  camera’s  electronic  chiller  to  5°  C 

Flat-Field  Correction 

—  Photograph  uniformly  illuminated 
surface 

—  Divide  images  by  flat-field  images  to 
correct  for: 

•  Vignetting  (up  to  60%  brightness 
variation  from  center  to  corner  in 
non-corrected  images) 

•  Non-uniform  pixel  response 

•  Dust  particles  on  CCD 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  unlimited. 


Response  (Counts) 


Relative  Radiance  Calibration 


•  Scale  images  through  relative 
radiance  calibration 

—  Integrating  sphere 

—  Quartz  tungsten  lamp  with  NIST 
traceable  irradiance  curve 

•  Calculate  and  Qg23  /  Q542 

•  Characterize  pixel  response  to  ensure 
exposures  are  within  the  linear  regime 


R  = 


=  c 


ens 


/ic^pixel  5'2 


At 


^1  _  ^1  Qx2  ^2  ^^2  fx2 
^2  ^2  Qxi  ^1  fxi 

C  =  combined  physical  and  geometrical  constants 
n  =  CCD  count  value 

=  fractional  transmittance  of  ND  filter  at  a  given  wavelength 
=  CCD  quantum  efficiency  at  a  given  wavelength 
A  =  wavelength 
h  =  Planck  constant 
c  =  speed  of  light 
/tpixei  =  area  of  CCD  pixel 
/\iens  =  area  of  lens 
s'  =  distance  between  lens  and  CCD 
At  =  exposure  time 
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Abel  Inversion  Technique 


Forward 

Transform 


Inverse 
T  ransform 


1  dF(y)  dy 

njy  dy  ^y2  _  j-2 


Inverse  Abel  Transform  is  not  practical 
to  solve  directly  with  numerical  means 

—  Singularity 

—  Requires  derivative  of 
experimental  data 

Use  “Fourier  Method”*  to  approximate 
f(r)  as  a  cosine  expansion 

Nu 

fir)  =  ^  Anfn  (r) 

Ni 

fair)  =  1, 

) 

Choose  frequency  limits  to  balance 
noise  filtering  and  detail  retention 

Least  squares  fit  of  forward  transform 
with  experimental  data  determines 
values  of 


fnir)  =  1-  (-l)"cos  (nn^ 


‘Pretzier,  G.,  Jger,  H.,  Neger,  T,  Philipp,  H.,  and  Woisetschlger,  J.,  “Comparison  of  Different  Methods  of  Abel  Inversion 
Using  Computer  Simulated  and  Experimental  Side-On  Data,"  Zeitschrift  fr  Naturforschung,  Vol.  47a,  1992,  pp.  955-970. 
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Operating  Condition 

V,{V) 

Power 

(W) 

(mg/s) 

m.  (mg/s) 

(mg/s) 

P  (Torr-Xe) 

P  (Pa-Xe) 

Nominal 

300 

4.61 

1383 

5.15 

0.40 

0 

2.28E-05 

3.04E-03 

High  Back.  Press. 

300 

4.57 

1371 

5.15 

0.40 

4.22 

4.55E-05 

6.07E-03 

Low  Dis.  Voltage 

210 

4.60 

966 

5.15 

0.40 

0 

2.09E-05 

2.79E-03 

High  Dis.  Voltage 

390 

4.56 

1778 

5.15 

0.40 

0 

2.08E-05 

2.78E-03 

Low  Prop.  Flow  Rate 

300 

2.90 

870 

3.60 

0.28 

0 

1.52E-05 

2.03E-03 

High  Prop.  Flow  Rate 

300 

5.99 

1797 

6.69 

0.51 

0 

2.74E-05 

3.66E-03 
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150mm  100  50  0  -50  -100  -150 


Comparison  of  Ion  and  Neutral  Emission 

Radiance  (arb.  unit) 


Z 

0  50  100  150  200mm 


o 

lO 


o 

o 


z 

0  50  100  150  200mm 


z 

100  150  200mm 


542  nm  (ion)  823  nm  (neutral)  Afi(%)  =  x 

^823:Wom 

100 

Neutral  emission  is  greater  than  ion  emission  throughout  plume,  especially  at  the 
cathode  (26x  greater  radiance) 

Highest  regions  of  relative  ion  emission  found  at  channel  exits  and  central  plume 
core,  lowest  regions  outside  of  plume  beam 
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150mm  100  50  0  -50  -100  -150 


Z 


542  nm  (Ion)  Emission  Double  BG  Press 

Radiance  (arb.  unit) 


Z 


Z 


50  100  150  200mm 


0  50  100  150  200mm 


Nominal  (P  =  2.3e-5  Torr)  High  Pressure  (P  =  4.6e-5  Torr)  A/?(%)  =  (^542:dp  R542-.Nom)  ^ 

^S42\Nom 

•  Higher  thrust  measured  as  background  pressure  increases  due  mostly  to  less  plume  divergence*  1 00 

•  Asymmetrical  region  of  high  emission  at  z  =  20,  y  =  75  mm 

•  Lower  emission  in  plume  core  and  near  exit  plane,  higher  elsewhere 

—  Maybe  less  beamlet  divergence  from  channel  exits,  from  upstream  shift  of  ion  acceleration  region 

•  Similar  emission  at  channel  exit. 


*Diannant,  K.  D.,  Liang,  R.,  and  Corey,  R.  L.,  “The  Effect  of  Background  Pressure  on  SPT-100  Hall  Thruster  Performance,"  Proceedings  of  the  50th  AIAA/ASME/SAE/ASEE  Joint  Propulsion 
Conference,  2014,  AIAA-201 4-3710. 
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150mm  100  50  0  -50  -100  -150 


Nominal 


Double  BG  Pressure 


A/?(%)  = 


(^823:DP  ^823:JVo?n) 

^823:JVom 


X 


100 

•  Similar,  but  more  intense  asymmetrical  region  of  high  emission  at  z  =  20, 
y  =  75  mm  as  ion  emission 


•  Broader  low  emission  region  in  plume  core  beyond  z  =  50  mm 
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Radiance,  F(y)  (arb.  unit) 


Radial  Profiles:  Double  BG  Press,  vs.  Nominal 
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Abel  Inversion  of  Radiance.  f(r)  {arb.  unit) 
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542  nm  (Ion)  Emission  Variations 
Radiance  (arb.  unit) 
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•  Greater  emission  in  plume  beyond  z  =  100  mm,  lower  upstream 


Maybe  beam  interaction  is  occurring  farther  downstream  with  the  faster  ions 
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Similar  changes  compared  to  ion  emission,  but  with  broader  plume  emission  features 
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Radiance,  F(y)  (arb.  unit) 


Radial  Profiles: 
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Abel  Inversion  of  Radiance,  f(r)  (arb.  unit) 
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542  nm  (Ion)  Emission  m  Variations 
Radiance  (arb.  unit) 
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Dramatic  increase  in  emission  for  high  propellant  flow  rate  (note  color  scale  change) 

—  Highest  increase  directly  downstream  of  center  magnetic  pole  and  outside  of  channel 

—  More  collisions  with  unionized  propellant 
Spherical  emission  shape  different  from  other  cases 

—  Scattering  collisions  diffusing  plume 
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823  nm  (Neutral)  Emission  m  Variations 

Radiance  (arb.  unit) 
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Similar  changes  to  plume  emission  shape  as  seen  for  542  emission 
—  High  increase  in  emission  directly  downstream  of  channel  exit 
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Conclusions  and  Future  Work 


*  Conclusions 

o  Optical  Technique 

•  High  degree  of  repeatability 

•  Exposures  of  separate  regions  of  plume 
scaled  well  for  composite  image  increasing 
dynamic  range  of  data 

•  Inverse  Abel  inversion  analysis  technique 
performed  well  overall 

—  Resolved  emission  detail  related  to 
acceleration  channel  geometry 

—  Similar  neutral  distribution  shape 
across  channel  as  seen  in  LIF  density 
measurements  SPT-140  near-field 
plume* 

—  A  few  problem  areas  near  r  =  0  where 
f(r)  <  0 

o  Plume  Analysis 

•  Distinct  changes  in  emission  structure  of 
plume  observed  for  various  operating 
conditions 


•  Future  Work 

□  Model  Comparison 

•  Perform  measurements  to  match  a 
numerical  model  of  the  near-field  plume 
that  includes  CRM 

•  Compare  experimental  and  numerical 
results  for  model  validation 

•  Examine  scaling  of  emission  with  density  in 
numerical  results  to  see  if  photographic 
density  estimates  are  appropriate 

□  Probe  Comparison 

•  Directly  measure  density  and/or 
temperature  with  a  probe  in  the  near-field 

•  Compare  results  with  emission  data  to  see 
if  density  can  be  estimated  from  emission 
where  temperature  gradients  are 
adequately  low 

•  Photograph  probe  in  plume  to  examine  the 
perturbation  as  function  of  emission 

□  Apply  experimental  technique  with  high 
speed  camera  to  study  breathing  mode 


*Crofton,  M.  et  al.,  “Neutral  Density  in  the  SPT-140  Near-Field  Plume,"  Proceedings  of  the  33rd  International  Electric 
Propulsion  Conference,  2013,  IEPC-201 3-399. 
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Forward  Abel  Inversion 


hiy)  =  2 


Forward  Abel  Inversion  becomes 
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/W  ^ - jd'' 
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H(y)  =  2  E  A,  \m—; 

«  =  J  |/j 

dr 


Inverse  Abel  Inversion 


/w-i  j 


Set  each  integral  term  numerically 

R 

r  dr 


r  my)  dy 

Kiy)  = 
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Lif) 

J  \/y^-r^ 
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} 

\ 

Approximate  f(r)  as  cosine  func.  expansion 

m=  I  AJ„{r) 

F|=iVj 

/oW  =  1  > 


/„(0  =  l-(-ircos|n;r- 


Perform  least  squares  fit  to  fit  amplitude  values  A„ 

n  =  jVj  \  lt=l  /  ic=l 

Vm:  Ni<m< 


N 


Jlr=  1 
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